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Circular  Array  of  Dipoles 
Above  a  Perfectly  Conducting  Cylinder 


1.  INTRODUCTION 

Recent  developments  in  low-inertia  phase  shifters  and  electronic  switches1 
have  aroused  interest  in  large,  circular  arrays  on  conducting  cylindrical  surfaces 
for  electronic  agile  beam  positioning  or  radar  resource  optimization,  invariant 
beam  characteristics  over  the  entire  360°  azimuthal  scan  sector,  and  frequency- 
independent  beam  pointing  direction.  A  low-sidelobe,  high-gain  aperture  design 
for  a  circular  array  of  dipoles  surrounding  a  perfectly  conducting  cylinder  requires 
a  knowledge  and  control  of  mutual  coupling  parameters  in  the  array  environment. 2 

3 

Carter  in  1943  outlined  a  rigorous  solution  of  Maxwells'  equations  for  a  dipole 
near  a  long  cylinder  (see  Figure  1),  and  tabulated  formulas  for  the  radiation  pat¬ 
terns  for  three  different  circular  array  configurations.  Carter  used  the  reciprocity 
theorem  and  the  k  town  solution  for  scattering  from  a  cylinder  to  obtain  far-zone 

(Received  for  publication  30  September  1980) 

1.  Provencher,  J.H.  (1910)  A  survey  of  circular  symmetric  arrays.  Phased 

Array  Antennas  (Oliver  and  Knittel,  Artech  House,  Inc.,  Dedham,  Mass.) 
page  292. 

2.  Hessel,  A.  (1970)  Mutual  coupling  effects  in  circular  arrays  on  cylindrical 

surfaces-«perture  design  implications  and  analysis.  Phased  Array  Antennas. 
(Oliver  and  Knittel,  Artech  House,  Inc.,  Dedham,  Mass. )  page  273. 

3.  Carter,  P. S.  (1943)  Antenna  arrays  around  cylinders.  Proceedings  of  IRE, 

Vol.  31,  December, 
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Figure  1.  Dipole  Near  a  Long  Cylinder 
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patterns.  Lucke  •  '  used  Green's  function  method  to  yield  expressions  for  the 
fields  in  terms  of  an  integration  in  the  complex  plane.  The  problem  was  also  solved 
by  Harrington  *  1  by  converting  a  three-dimensional  radiation  problem  having  cy¬ 
lindrical  boundaries  into  a  two-dimensional  problem  by  applying  a  Fourier  trans- 
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form  with  respect  to  the  cylinder  axis,  Z.  Noif  also  used  the  Fourier  transform 
method,  and  by  means  of  a  limiting  process,  obtained  an  approximate  form  for  a 
dipole  near  a  finite  cylinder. 

In  this  report,  the  problem  of  finding  the  electromagnetic  field  and  admittance 
of  a  circular  array  of  dipoles  near  and  parallel  to  an  infinite  conducting  cylinder 


4.  Lucke,  W.  (1949)  Electric  Dipoles  in  the  Presence  of  Elliptical  and  Circular 

Cylinders,  Keport  No.  1,  Project  188,  Stanford  Research  Institute, 

Stanford,  Calif. 

5.  Lucke,  W.  (1951)  Electric  dipoles  in  the  presence  of  elliptical  and  circular 

cylinders,  journal  of  Applied  Physics,  Vol.  22,  No.  1,  January. 

6.  Harrington,  H.  F.  and  LePage,  W.  R.  (1949)  A  Study  of  Directional  Antennas  fern 

DIF  Purposes,  Report  1,  Department  of  Electrical  Engineering,  Syracuse 
University,  Syracuse,  NY,  September. 

7.  Harrington,  K.  '.  (19(>1)  Time -Harmonic  Electromagnetic  Fields,  McGraw-Hill, 

New  York,  NY,  Chapter  3. 

8.  Neff,  H,  P. ,  Hickman,  C.  E. ,  and  Tillman,  J.D.  (1964)  Circular  Arrays  Around 

Cylinders.  Report  No.  7,  Department  of  Electrical  Engineering,  University  of 
Tennessee,  Contract  AF19(628)-288,  June. 
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has  been  mathematically  formulated  Into  Integral  equation*  via  the  Fourier  trans- 
form  method,  which  relates  the  illumination  function's  voltages  to  the  resultant 
current  distribution  due  to  the  mutual  coupling.  We  use  the  computer  to  numerically 
solve  the  m  nonhomogeneoua  integral  equations,  using  numerical  approximations 
for  the  integrals  and  complex  matrix  techniques  to  solve  the  simultaneous  equations. 
Such  an  analysis  should  account  for  both  the  space  wave  with  Ita  grating  lobe  and 
the  creeping  wave  resonances  on  the  curved  surface,  and  have  an  impact  on  the 
development  of  optimum  method!  for  determining  techniques  for  radiation  pattern 
generation  with  very  low  sidelobes,  optimum  spacing  and  a  minimum  number  of 
array  elements  (see  Figure  2).  We  snail  investigate  the  effects  of  mutual  coupling 
for  a  circular  array  located  a  quarter  of  a  wavelength  above  a  perfectly  conducting 
cylinder  having  a  radius  of  two  wavelengths.  This  simple  case  was  chosen  to 
easily  investigate  the  m  simultaneous  Integral  equations  and  to  lay  the  foundation 
for  future  investigations  into  the  more  practical  case  of  a  perfectly  conducting 
cylinder  with  radius  of  approximately  10  wavelengths. 


objective:  TO  OBTAIN  360-DEGREE  scan  capability  for 
LOW  SIDELOBE  ARRAYS  OF  ELEMENTS 


1.  MUTUAL  COUPLING  EFFECTS 

2.  OPTIMUM  ILLUMINATION 
FUNCTION 


I  LOW  SIOELOBE  AZIMUTHAL 
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2.STACKING  OF  CIRCULAR 
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APPLICATIONS.  UNATTENDED  RADAR,  GROUND  BASED  TACTICAL  RADAR 


Figure  2.  Circular  Arrays  Around  Cylinder,  Optimization  Approach 


2.  MATHEMATICAL  ANALYSIS 


The  circular  antenna  array  to  be  considered  (see  Figure  3)  consists  of  Identi¬ 
cal,  parallel,  cylindrical  dipoles  equally  spaced  around  the  circumference  of  a 
circle,  which  is  concentric  to  the  perfectly  conducting  cylinder  which  the  array  of 
dipoles  surrounds.  Only  center- fed  dipoles  with  a  half-length  X/4  and  a  radius  b 
will  be  used  as  elements.  The  axial  center  of  each  dipole  is  perpendicular  to  the 
plane  of  the  circle.  Each  of  the  elements  of  the  array  is  thus  in  the  same  geomet¬ 
rical  environment.  The  radius  of  the  perfectly  conducting  cylinder  is  a,  and  the 
radius  of  the  circular  array  of  dipoles  is  rQ  (see  Figure  4).  The  expression  tor 
the  vector  potential  A  (r,  4,  z)  has  been  developed  for  a  circular  array  of  m 
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dipoles  located  about  a  perfectly  conducting  cylinder. 

For  a  current  density  j  located  above  a  perfectly  conducting  cylinder  (see 
Figure  5), 

V2  X  +  k2  J  =  -„o7  .  (1) 

For  a  unit  dipole  at  (r  .  4  ,  z  ), 
o  o  o 

-  .  *<r  -  r  )6(4  -  4  >6(z  -  *  > 

j  =  z  - 2 - _2 - 2_  .  (2) 

Therefore,  we  need  only  be  concerned  with  the  Az  component: 

00 

f(z)  =  37  J  F(o.  r.  4)  eloZ  do  (3) 

where 

00 

F(o,  r,4)  =  J]  am(ff.  r)eim*  .  (4) 

I’jl  =  -00 


9.  Fante,  Dr.  Ronald  L.  (private  communication). 
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CYLINDRICAL  BOUNDARY 
CONTAINING  CURRENT 
ELEMENTS 


Figure  4,  Perspective  View  of  Dipole,  Conductive 
Cylinder  Geometry 
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sun 


Figure  5.  Geometry  for  the  Development  of  Integral  Equations  for  a  Circular 
Array  Around  a  Perfectly  Conducting  Cylinder 


From  Eqs.  (3)  and  (4)  we  get 


In  cylindrical  coordinates,  the  laplacian  is 


.13 

V  A'-r  5F  ^ar 


Substitute  Eqs,  (G),  (5),  and  (2)  into  Eq,  (1)  to  yield 


7*  £ 


,i«z  J  1  3 


r  5r  I  9r 


/.2  2  m"  \ 

+  k  -  a - k—  \  a 


6(r  -  r  )6(*  -  6  )*<z  - 


T 
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'  i  ijjjpuv  jnpRntRP"  ■ 


B  -  D 
n  n 


c  j  (x) 

-ft^r  • 


Substitute  Eq.  (IS)  into  Eq.  (14)  to  yield 


Define  A  =  J_(x>  H’(y)  -  JMx)  H(x)  (Wronakian),  then 
n  n  n  n 


-M  H(X) 
o  n 

2*  X  A 


Also  from  Eq.  (15), 


(15) 


(10) 


-n  J  (x) 

n  .  r>  -  °  n 

n  un  2*X  A 


(17) 


Now  we  need  to  satisfy  the  boundary  conditions  at  the  surface  r  -  a  of  the  cylinder, 
which  is  that  the  electric  field  tangential  on  its  cylindrical  surface  must  be  zero. 


E  - 


Jw<o% 


VX  vx  A 


.  /  9A  \  .  9 A 

**A  *  fF  ar  -♦  FT2- 


VX  VX A  =  r 


92A  .  ,  92A 

_ £  ,  T  1  Z  ,  * 

drdz  *  r  5zdj> 


1  9 

(r^ 

.  1 

32A 

z 

"r  Jr 

V  / 

T 

3(j>^ 

r 

therefore 


dzdtj) 


=  0 


r=a 


(18) 
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14 


d  ■  'c",»'y- 

Dn  ’  ~H^y» 


Upon  substituting  Eq.  (16)  into  Eq.  (24)  and  rearranging  terms  we  get 


,  »oHn(x»  (>n'»\ 

V  ■  Tors'  j 


Using  Eq.  (25)  in  Eq.  (17)  we  get 


Since  *  X  A  =  2i  we  get 


B„  -  jh^{V’‘'h»w-  W'} 

_  ..  K«r,M 
Cn  "  - 3 


D 


n 


H  (x) 
o  n 

4i 


Jn<y> 

Hn(y) 


(24) 


(25) 


(26) 


(27) 


(28) 


(29) 


where 


rQ  =  radius  at  which  dipole  is  located 
a  =  radius  of  cylinder. 
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fo*1  rsr^fi 


e  T  *.w4  ]  '*'''«*■  ±  «“<vb/r”‘v 

m*l  -L/2  ^  ■»  n*-« 


*  Bn<a)Hn^2  ■  °2  r)  *  T  [ck  cot  **«  +  ^L-tn  <*M  ]  (33> 

where  the  K**1  dipole  is  at  (r,  i,  a).  Thus,  for  an  m -dipole  circular  array,  we  have 
in  simultaneous  equations  to  solve.  Let  a  =  */20,  then  Eq.  (33)  becomea 


m  L/2  go 

AE  /  d*°  w  / 

ml  -L/2  -« 


eio(t/23-ao)  ^  £  ein4k+b/ro-dm) 
n=-*e 


X  12  -  o2  r)  = 


n  n' 


(34) 


*z 


Substitute  Eq.  (27)  into  Eq.  (34)  and  let  I  (a  )  =  A  cos  then, 

mom  L 


m 


L/2 

i?E*.  / 

m^l  -L/2 


dz  cos 

o 


az  ® 
*o  f 

L  J 


dw  e 


1  £  ^n«k.b/,o-4m) 


n  -go 


rj  (*)II  (y).H  Ix)J(y)  /- 5 - 5  T  -4»V.. 

n/  I  -13 - 0 _ n _ Q__  ..  <*/ 1.2  2  _v  I  k 

L  H  (y)  Hn(Vk  -*  r)J  -  -Te¬ 


rn) 


Since  all  dipoles  are  equidistant  from  the  perfectly  conducting  cylinder,  r  *  r 


then, 


Hf'/  -  o2  r)  I  I  (S  k“  -  J  r  )  =  H 


o  n 


(x) 
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which  from  symmetry  becomes 


m 


E 

m  =  l 


eln(*k +b/ro"*m>j? 

0 


cos 


4 


X 


J_(x)  H  (x) 
n  n 


hW£) 

Hn«» 


iXiii 

*M„C 


(3?) 


2.1  Saif  mi  Mutual  Adauttaec* 

Now,  let 


« 

Tkm(“>  =  5Z  e 

n=0 


in<Vb/ro*V  J  da  -- 


Lo 

cos  X_  coa 

o  4 


X 


J(x)H  <x) 
n  n 


Hj(x)Jn<y) 

H„W 


then 


m 

E 

m-i 


AmTkm(o> 


'JjsL 


Now  we  shall  consider  the  M  integral  equations  by  taking  a  look  at  the  equations 
yielded  for  each  base  voltage  (K  -  l,  ...  m)  simultaneously.  We  get 


A1T11  +A2T12  +  AmTlm  ’  JTc  VI 
AlT21  .  A2T2J  .  ...  AmT2m  ■  V2 


A,Tml  +  A,Tm,  +  ...  A  T  -  -k*  Vm 
I  ml  2  m2  m  mm  u  L  m 


(38) 
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In  a  matrix  notation,  Eq.  (38)  ia  equal  to 


Ttl  TU 


Tml  Tm2 


* 

*  - 

P 

Tlm 

Ai 

*« 

T2m 

A2 

V, 

• 

• 

s 

* 

• 

T 

mm 

• 

A 

m 

m 

where 


(38) 


Therefore, 

f  X  =  V  . 


Given  values  for  V^.  V2,  ....  V  and  having  numerically  approximated  the  inte¬ 
grals  Tkm  (m  =  1 . mj  K  =  l . m),  we  compute  the  values  for  the  Am's 

by  factoring  the  matrix  T  into  the  L  -  U  decomposition  of  a  row -wise  permutation 

of  T  and  solving  the  systems  of  Eq.  (37).  Once  the  values  for  the  A  's  have  been 

m 

found,  we  can  then  determine  the  mutual  and  seif  admittances  of  the  circular  array. 
Equation  (33)  gives  the  relationship  between  the  currents  and  voltages  or  the  array 
elements.  For  an  m-element  array. 


YUV1+Y12V2  +  -  +  YlkVk  +  -  "YlmVm=1l 
Y21Vl+Y22V2  +  "‘  +  Y2kVk+~*  +Y2mVm  =  l2 


(40) 


YmVl  +  Ym2V2 


+  YmkVk  + 


+  V  =  I 
mm  m  m 


where  Y^,  Y22.  •••  Ykk,  ...  Ymm  are  the  self  admittance  of  the  respective  ele¬ 
ment,  Yflb  is  the  mutual  admittance  between  a  and  b,  is  the  applied  voltage  on 
the  K1*1  element,  and  1^  is  the  current  in  the  element. 

For  an  array  with  a  single  active  element,  which  for  simplicity's  sake  we  let 
be  element  no.  1,  every  VK  where  K  =  1  is  equal  to  zero.  Thus,  Eq.  (16)  becomes 
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*1  *  Yllvl 
*a  *  Y2ivi 
*3  -  Y31V1 


(41) 


We  let  Vj  «  1  end  ta  the  magnitude  and  phase  of  the  current  on  element  K, 
Therefore, 

A1  =YU 
A2  =  Y21 


Am  *  Yml 


(42) 


Thus,  the  coefficient  Aj  la  the  self  admittance  of  the  active  element  no.  1  and  the 
coefficient  AVl  where  K  Is  unequal  to  one,  is  the  mutual  admittance  between  the 
first  and  K  element. 


2.2  Numerical  Reaulta 

The  simple  case  of  a  =  2A.  and  r  =  a  +  A./4  (see  Figure  7)  was  investigated 
and  the  mutual  and  self  admittances  were  determined  through  use  of  computer 
programs  to  numerically  solve  the  m  nonhomogeneous  simultaneous  equations 
(Eq.  (37)).  Equation  (37)  becomes 


m 

oc 

E 

Am  E 

m  =  1 

n=0 

in(dk+b/r0’4m) 


BESS(n)  = 


“VKL 

2pqc 


(43) 
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where 


k  » 

BESS(n)  *  J  doQ(o)  ReFS(n,  a)  +  ^  do  Q(o>  ReFB<n,  o) 
o  k+e 


♦ 


do  Q(o)  lmFS(n,  o)  + 


« 

; 


do  Q(o)  ImFBin,  o) 


} 


(44) 


Q(o)  « 


(45) 


ReFS(n.o)  «  JR(x)Jn(x) 


Jn(x>Jn(y>  -  Yn(x>Jn(y>  +  aJB(«)Vn(x)Jn<y)Yn<y) 


Jj(y>  +  Y^(y) 


(46) 


ReFB(n.o)  *  IR(x)KR(x) 


(47) 


2Jn(x)Yn(x)j;(y)  -  JR(x)JR(y)YR(y)  4  Y*(x)JH(y)Yn(y) 


ImFS(n.a)  •  J(x)Y  <x) - s - 

n  n  Jn<y>+Yn<y) 


(48) 


2  J  K?(x)In(y) 1 
ImFB(n.a)  ■  |  \~\^r 


(49) 


For  Eq.  (43)  to  be  solvable,  the  expression 


'A  in(*i+b/r-4) 

>  f  e  K  0  m  BESS(n) 

n=0 


(50) 


must  converge  to  aero  for  large  n,  otherwise,  the  appropriate  coefficients  for  each 
Am  (m  *  1,  2,  . . .  m)  cannot  be  determined.  Let 


SR(n)  *  /Qk)ReFS(n,o)  do 
SI(n)  =  /Q(o)ImFS(n,o)  tfe 
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BR(n)  *  /Q(o)ReFE(n,er)  do 
BI(n)  -  / Q(or)ImFB(n, or)  db 

The  coi  trlbutions  to  BESS(n)**  of  SR  and  BI  rapidly  diminish  for  large  order  n. 

BR  and  SI  decrease  very  slowly  in  absolute  vdtue  for  lrrge  order  n.  Figures  8,  9, 
10,  and  11  show  the  values  of  SR,  SI,  BR  and  BI  respectively  for  order  n  =  0  to 
700.  BR  and  SI  approach  zero  very  slowly  for  large  order  n,  while  SR  and  SI  con¬ 
verge  rapidly  to  zero  for  large  n.  Figure  12  shows  the  value  BST.1M(n)  where 

BSUM(n)  *  J (SR(n)  +  BR(n))2  +  (SJ(n)  +  BI(n))2  *  |BESS(n)j 


and  ^  *  0°  (worst  case).  BSUM  is  shown  to  slowly  converge  to  zero  for 

large  order  n.  Due  to  time  limitation  in  computer  usage,  n  -  0  to  700  was  chosen; 
however,  for  more  reliable  results,  n  *  0  to  2000  may  be  required,  especially  for 
the  worst  case  *  0°. 

Once  l.\c'  BESS(n)  were  determined  for  n  =  0  to  700,  these  values  were  used  in 

Eq.  (50)  to  determine  the  coefficients  for  the  A_  's  in  the  set  of  simultaneous  non- 

m 

homogeneous  equations  (Eq.  (37)).  Each  resultant  matrix  was  found  to  be  of  the 
form 


T1  Tm  *  *  *  T2 

iV 

V1 

T2  T1  •••  ***3 

A2 

V2 

•  i  •  • 

A3 

V3 

•  «  •  • 

A4 

V4 

•  •  •  • 

• 

• 

T  T  t  •  T. 

A 

V 

m  m-1  1 

m 

m 

(51) 


Thus,  due  to  symmetry,  Eq.  (37)  need  only  be  evaluated  once  for  any  m  =  i,  2,  ... 

m  to  determine  all  coefficients  for  the  A  's  in  the  equation. 

m 

Results  of  the  self  and  mutual  admittances  for  a  circular  array  of  cylindrical 
dipoles  a  quarter  of  a  wavelength  above  a  perfectly  conducting  cylinder  of  radius 
r  -  2A  were  calculated.  See  Figure  12  for  array  sizes  of  m  =  2  to  22.  Figure  13 
shows  the  self  admittance  for  the  active  element,  Figure  14  shows  the  mutual  ad¬ 
mittance  between  dipole  no.  1  and  dipole  no.  m.  Due  to  symmetry,  these  two 
should  be  the  same.  The  inconsistency  may  be  the  result  of  not  taking  large  enough 


11.  Amos,  V.  E, ,  and  Daniel,  S,  L.  (1977)  AMOSLIB,  A  Special  Library  Version. 
Sandla  Laboratories  (SAND  77-1390). 
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order  n  in  Eq.  (50)  to  have  BESS(n)  approximately  equal  to  zero.  Figure  15  shows 
the  mutual  admittance  between  dipole  no,  1  and  the  adjacent  dipole  no.  m.  Fig¬ 
ure  16  shows  the  relative  currents  for  the  dipoles  for  arrays  of  size  m  =  2,  4,  7, 
10,  13,  16.  Once  again,  for  m  a  12,  asymmetries  occur  which  may  be  due  to  not 
taking  high  enough  order  n  for  the  evaluation  of  BESS(n)  in  Eq.  (50)  for  the  deter¬ 
mination  of  the  coefficients  for  the  Am's. 


Figure  7.  Array  of  Cylindrical  Dipoles  About  a  Perfectly  Conducting 
Cylinder,  Simple  Case,  a  =  2X,  p  =  rQ  =  a  +  X./4 
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Figure  10.  BR(n)  =  f  Q(a)  ReFB(n,  o)  for  n  =  0  to  700 


Figure  11.  BI(n)  -  /Q(o)  ImFB(n,  a)  for  n  =  0  to  700 
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Figure  16.  Relative  Currents  For  the  Circular  Array  of 
Dipoles  Above  a  Perfectly  Conducting  Cylinder  For  Active 
Element  Dipole  No.  1,  m  =  2,  4,  7,  10,  13,  16 
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Figure  16.  Relative  Currents  For  the  Circular  Array  of 
Dipoles  Above  a  Perfectly  Conducting  Cylinder  For  Active 
Element  Dtpole  No.  1,  m  =  2,  4,  7,  10,  13,  16  (Cont. ) 
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S.  CONCLUSION 


It  has  been  shown  that  numerical  techniquea  for  approximation  of  the  aolution 
of  the  m  nonhomogeneous  aimultaneous  equations  (Eq.  (37))  provide  useful  and  ana¬ 
lytical  data  in  the  determination  of  characteristics  of  a  circular  array  above  a 
perfectly  conducting  cylinder.  From  Eqs.  (37)  and  (SO),  it  is  seen  that  the  value 
BESS(n)  must  converge  to  zero  for  large  order  n  to  determine  the  coefficients  of 
the  Am'a  in  Eq.  (39).  Due  to  the  symmetry  of  the  coefficients,  as  seen  in  Eq.  (51), 
the  values  Tj  j,  T^.  •  •  •  Tm  need  only  be  calculated  once  for  a  given  array  of  slxe 
m.  Also,  from  the  expression  for  BESS(n),  Eq.  (50),  BESS(n)  is  independent  of  m, 
the  number  of  dipole  elements  In  the  circular  array  being  considered.  BESS(n)  is, 
however,  dependent  on  the  values  of  the  radius  of  the  cylinder  and  the  frequency 
(Wavelength).  Consequently,  tables  can  be  made  of  the  BESS(n)  values  for  n  =  0 
to  N,  where  N  is  of  order  n  large  enough  to  force  BESS(n)  to  zero.  These  tables 
would  enable  calculation  of  the  mutual  and  self  admittances  for  an  array  of  arbitrary 
size  m  about  a  perfectly  conducting  cylinder  of  radius  r'.  An  understanding  of  the 
effects  of  mutual  coupling,  taking  into  consideration  both  the  space  wave  and  the 
creeping  wave,  is  necessary  in  the  synthesizing  of  an  aperture  distribution  to  give 
the  best  fit  to  a  specified  radiation  pattern,  and  for  the  determination  of  optimiza¬ 
tion  techniques  for  circular  array  antennas. 
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